Intense terahertz laser fields on a quantum dot with Rashba spin-orbit coupling 
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We investigate the effects of the intense terahertz laser field and the spin-orbit coupling on single 
electron spin in a quantum dot. The laser field and the spin-orbit coupling can strongly affect 
the electron density of states and can excite a magnetic moment. The direction of the magnetic 
moment depends on the symmetries of the system, and its amplitude can be tuned by the strength 
and frequency of the laser field as well as the spin-orbit coupling. 
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The study of semiconductor based opto-electronics, 
which combines the ultrafast electronics with the low- 
power optics, has been fruitfully carried out for years. 
Schemes based on semiconductor quantum dots (QD's) 
have been proposed and actively studied in order to real- 
ize the more ambitious goal of quantum opto-electronic 
devices, which provide interface between quantum bits 
(qubits) and single-photon quantum optics and can per- 
form quantum information processing and communica- 
tion on one chip. Single spins in QD's are natural 
candidates^'^'^ for qubits as they have long coherence 
times and can be manipulated both electronically and op- 
tically. There are many theoretical works on spin relax- 
ation/decoherence and manipulation in QD's.'''^'^'^ Ex- 
perimental realization of the electrical/optical generation 
and read-out of single spin in QD have also been carried 
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The spin-orbit coupling (SOC) plays an essential role 
in the dynamics of electron spins in QD's. It enables 
the electrical manipulation by tuning the gate voltage 
but also leads to the spin relaxation/decoherence to- 
gether with the different scattering mechanism such as 
the electron-phonon scattering.^ In addition to the above 
effects, the SOC in QD's causes a spin splitting of a few 
meV, the order of terahertz (THz), and therefore is ex- 
pected to affect the response of electrons in QD's to the 
THz laser. The effect of the SOC on far-infrared op- 
tical absorption spectrum has been discussed in detail 
by Rodriguez et al}^ Until recently, the study on the re- 
sponse to the THz radiation has been focused on the spin- 
unrelated problems such as the dynamic Franz-Keldysh 
effect (DFKE), "'12.13.14 the AC Stark e^eci^^'^^ and 
the photo-induced side-band effect. ^^-^'^'i^'i^'^'''^! Very 
recently Cheng and Wu discussed the effects of intense 
THz laser field on a two dimensional electron gas (2DEG) 
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with the Rashba SOC.^^ It is demonstrated that the laser 
field can significantly modify the density of states (DOS) 
of the 2DEG and induce a finite off-diagonal density of 
spin polarization which indicates the strong correlation 
of different spin branches. Later Wang et al. discussed 
the time-dependent spin-Hall current in the 2DEG driven 
by an intense THz field. ^"^ How the intense THz laser af- 
fects the QD's in the presence of the SOC needs to be 
further revealed. In this paper, we study the effect of in- 
tense THz laser on semiconductor QD's with the Rashba 
SOC.24 

We consider a quantum dot grown in an InAs quantum 
well with growth direction along the z-axis. A uniform 
radiation field (RF) E(i) = E cos(i7i) is applied along the 
x-axis with the period To = if • By using the Coulomb 
gauge, the Hamiltonian is written as^^ 

H{t)^^ + V,{v)+Hso{^). (1) 

Here the momentum P — —i'V + eA{t) with A{t) = 
—Esm{nt)/n being the vector potential; m* is the elec- 
tron effective mass. The confining potential of the quan- 
tum dot is taken to be Vc(r) = ^'m*ujQr'^ in the x-y plane. 
An infinite-well-depth assumption is made along the z- 
axis and the well width is assumed to be small enough 
so that only the lowest subband is relevant. Hso is the 
SOC which is composed of the Rashba term^"* and the 
Dresselhaus term.^^ For InAs quantum well, the Rashba 
term is the dominant one. Hso(P,t) — aji[axPy — (jyPx\ 
with (T denoting the Pauli matrix and representing 
the Rashba SOC parameter which can be tuned by gate 
voltage up to the order of 4x10^^ eV-cm.^^'^* 

By employing the Floquet theorem, the solution to the 
Schrodinger equation with time-dependent Hamiltonian 
(1) can be written as,^^ 

^y^{Y,t) = e-»[£;™t-7sin(2nt)] 
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Here Egm = 4rn*'n'^ energy induced by the RF due 

to the DFKE;i2.i3 ^ = {'/>a(r)} can be any 

set of complete basis. {F^^} are the eigenvectors of the 
equations 



^ 2_^{an\Hj.\(3m)F^0 
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in which (r,t| an) 
and (an|7ijr|/3m) 



E </)a(r)e™"*, = H{t) - idt 

K?"" + mVL5c05nm with //" = 
7jr J^"dte^"'^^* H{t) representing the n-th Fourier com- 
ponent of the Hamiltonian. The eigen-values of the 
equations are solved in the region (— 0/2, f2/2].^^ In the 
case of the Rashba SOC, the non-zero terms are 



2m* 2 

{H-y = 



* 2 2 

+ -m Ur.r 



ieE 
2mQ 



Px - 



iaReE ^ 



20 



(4) 
(5) 



with p = — iV. 

We choose the complete basis {<^a} to be the eigen- 
states of the single spin in a QD without the RF and 

the SOC, whose Hamiltonian is described by the first two 
terms of Eq. (4). These states are identified as {| A^, / , a)} 
with 



{v\N,l,a) = 



x(r/a)l'le-(''/«)'/2Ll^l ((r/a)2)e''^Xa • (6) 



Here r = (rcos^, rsin0); a = y^hJmFuQ] lIJJ, is the 
Laguerre polynomial; rir = {N — \l\)/2\ and x<t is the 
eigenspinor of c^. With these basis functions one is able 
to write out the matrix Tij: and numerically diagonalize 
it.^° In this way, one obtains the eigenvalues e\ and the 
eigenvectors {F^^}. 

Without the RF and the SOC, the cigenstates of a sin- 
gle spin in the QD \N,l,(j) are 2 x (iV-|- l)-fold degenerate 
in angular momentum (index I = —N,—N + 2,--- ,N — 
2,N) and spin momentum (index cr =t, i), with energy 
being (A'' + l)a;o- The degeneracy is partially lifted by the 



SOC into {N + 1) states, each with 2-fold Kramers de- 
generacy. Although these states arc no longer cigenstates 
of az, they can be distinguished by the corresponding 
majority spin components as quasi-spin-up and -down 
states. The schematic of the lowest 12 states together 
with the possible transitions among them is shown in Fig. 
1. For both the small SOC and the small RF, the tran- 
sitions are mainly between the states in the same quasi- 
spin branch (the solid arrows in Fig. 1), since the RF 
does not flip the spin. At high RF and/or high SOC, the 
direct transitions (the dotted arrows) and multi-photon 
processes between different quasi-spin branches are im- 
portant, therefore the two quasi-spin branches become 
strongly mixed. The mixing can be roughly quantified 
by the quantity A = aaeE / {^l"^) which comes from the 
photon-assisted spin flip (the second term of H^) divided 
by the THz frequency.^^ The stronger the correlation be- 
tween different spin branches is, the larger A becomes. 
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FIG. 1: Schematic of the lowest 12 states of zero field Hamil- 
tonian H°, with all possible transitions. The transitions are 
indicated by solid arrows for quasi-spin conserving transitions 
and dotted arrows for quasi-spin flip transitions. 

The effects of the RF and SOC can be qualified by 
the density of states (DOS), pa.a, and density of spin po- 
larization (DOSP), pa-a-^^ These densities can be cal- 
culated from spectral functions which can be extracted 
from the Green functions: 



oo 

— OO 

OO oo 

= ^ / dte^-' j jdv ^ (r, hWx^ * (r, h) 

— OO — OO 

= J2 9{<Tu(r2;a,(3)F^^F^^*Jp{2'ycos{2flT))exp[i{n-m)QT] 

p,n,m.X,a.f3 

x6{u>-[ex+Eem-{2p + n + m)n/2]), (7) 



3 



where C and G° are retarded and advanced Green's 
functions respectively, T = (^i +t2)/2, t = ti — t2, cx and 
/? stand for {TV, I,(t}, Jp is the Bessel function of p-th 
order and 
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{N,l,a},l3 = {N',l',cT'}) 



(8) 



It is seen from Eq. (7) that these densities are periodic 
functions of T with period Tq. Symmetry analysis (time 
reversal symmetry and parity symmetry) further suggests 
that there is no spin polarization along the x- and z-axis. 
The only non-vanishing component of the spin polariza- 
tion is along the y-axis and reduces to zero if the RF is 
turned off. The spin polarization is an odd function of 
time p^ _^(T,uj) = — p„ _„{—T,uj), and therefore there 
is no overall spin polarization once averaged over time. 
On the other hand, the DOS is an even function of time 
and p|,|(T,w) = pi^i{T,u). Thus the RF and the SOC 
can affect the overall DOS, but cannot induce any spin 
polarization along the z-axis. In order to account for the 
effect of the RF and the SOC on the system, we use the 
average of p^^{T,uj) over a period to quantify the DOS 
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By numerically diagonalizing the matrix Tijr, one is 
able to obtain the coefhcients {F^^} and then calculate 
the DOS and the DOSP through Eq. (7). One can fur- 
ther obtain the time-averaged DOS by using Eq. (9) and 
calculate the magnetic moment along the y-axis My via 
proper evaluation of (py)?"^ To converge the DOS in the 
energy range between to 2.4wo, one has to use as many 
as 132 states from the lowest 11 major shells (TV = 0-10) 
since the energy levels of the QD are almost equidistant 
and hence lead to significant resonant overlap. In the 
following we present the numerical results of the time- 
averaged DOS and the magnetic moment My under dif- 
ferent conditions. In the calculation, the effective mass 
TO* = 0.0239mo with toq representing the free electron 
mass and the Rashba coefficient a/j = 3 x 10~^ eV • 
cm.^^^^^ The confinement potential ujq is chosen to be 5 

THz, which corresponds to a QD with diameter about 55 
nm.4'30,31 

In Fig. 2 we compare the time-averaged DOS with and 
without the SOC or the RF. The red (solid) curves are 
the DOS with both the SOC and the RF when E = 0.6 
kV/cm and = 4 THz. The blue (dotted) curve in Fig. 
2(a) is the DOS with only the SOC while the one in (b) is 
the DOS with only the RF. It is noted that in the figure 
the (5- function 8{x) in the DOS is replaced by the Gaus- 
sian function exp(— a;'^/2(T^)/(crA/27r) with a = 0.005^0- 
The DOS without the SOC and the RF peaks at = 
{N + l)ujo with value (iV + l)/{aV2^){N = 0,1,2, •••) 
(with the lowest two levels labelled by dashed lines in the 
figure). Each level has (A^-f l)-fold degeneracy as pointed 
out above. Nevertheless, by including only the SOC, the 
2-fold degeneracy of iV = 1 major shell is lifted and the 
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FIG. 2: (Color online) (a) The time-averaged DOS of QD 
under THz field with E = 0.6 kV/cm, = 4 THz (red solid 
curve) and _E = (blue dotted curve) for q_h = 3 x 10~^ 
eV-cm. (b) The time-averaged DOS of QD with Rashba SOC 
Qfjf, = 3 X 10~^ eV-cm (red solid curve) and an = (blue 
dotted curve) for E = 0.6 kV/cm, ^ = 4 THz. The black 
dashed lines correspond to the positions of and 2uo. 



ground state energy (iV = 0) is lowered by O.Icjq, as in- 
dicated in Fig. 2(a). Moreover, including only the RF 
gives rise to many peaks at the integer multiplications of 
ri, as can be seen in Eq. (7). It is seen from Fig. 2(b) 
that these peaks are equidistant. The peaks are gener- 
ally higher than those in the case with both the RF and 
the SOC, as the SOC causes additional splitting of the 
degenerate states. 

In order to further reveal the effect of the SOC and the 
RF on the system, we study the DOS under different RF 
and SOC. In Fig. 3, the DOS with different RF strengths 
(a) and Rashba coefficients (b) are plotted against en- 
ergy. It is seen from Fig. 3(a) that the first peak of the 
DOS gets a red shift from 0.75wo when E = 0.5 kV/cm 
to 0.68^0 when E = 0.6 kV/cm and then comes back to 
0.74a;o when E is further increased to 0.7 kV/cm. This 
is a clear demonstration of the combined effect of the 
AC Stark effect and the DFKE.^^ The AC Stark effect, 
which states that the energy difference becomes larger 
(smaller) if two energy levels are driven by an AC field 
with frequency below (above) the resonance. ^"'■^^ On the 
other hand, the DFKE states that the DOS has a blue- 
shift due to Eem as indicated in Eq. (7).^^^^^ In the cur- 
rent case, the frequency of the THz irradiation is 4 THz, 
which is smaller than the resonant frequencies of the tran- 
sitions from iV = to iV = 1 major shells even when the 
SOC is included. Therefore there is an AC Stark effect 
which shift the iV = major shell to red-side. This ef- 
fect is responsible for the red shifts as the DOS in the 
low energy range mainly consists of TV = major shell. 
Although both the DFKE blue-shift E^m and the AC 
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FIG. 3: (Color online) The time-averaged DOS (a) for E = 
0.5 (green dot-dashed curve), 0.6 (blue dotted curve) and 0.7 
(red solid curve) kV/cm, with an = 3 x 10~^ eV-cm and 
57 = 4 THz and (b) for ur = 1 (green dot-dashed curve), 
2 (blue dotted curve) and 3 (red solid curve) xlO~^ eV-cm, 
with E = 0.6 kV/cm and SI = 4 THz. 



Stark red-shift increase with the RF strength, the AC 
Stark effect is more important at low RF strength, and 
it saturates at high RF strength where the blue-shift due 
to the DFKE dominates. Moreover, at low RF inten- 
sity the peaks locate at e\ + E^rn — nQn/2 of a certain 
uq or replicas Sx + — [no + n)fl/2 with small n's. 
When the RF increases, the multi-photon processes be- 
come more and more important. As a result, the DOS 
becomes smoother as more and more multi-photon repli- 
cas appear. The effect of the SOC is plotted in Fig. 3(b), 
where the DOS under different Rashba coefffcient aj^ is 
plotted when the RF E = 0.6 kV/cm and 17 = 4 THz. 
It is seen that the first peak shifts to the red side as au 
increases. This is understood as the SOC contributes to 
the AC Stark effect. From Eq. (5), one notices that the 
AC Stark effect enhances as the SOC is increased. Differ- 
ing from the situation shown in Fig. 3(a), the blue shift 
Eem from the DFKE does not change since the RF is 
fixed here. Therefore the shift of the first DOS peak is 
monotonic. Moreover, as an increases the multi-photon 
processes between different spin branches become more 
important and the DOS becomes smoother. 

As said before, the DOSP is not zero when both the 
SOC and the RF are present. Due to the symmetry of 
the system, the only remaining spin polarization is along 
the y-axis. Consequently, the induced average magnetic 
moment can be calculated through the equation^^ 



M{T)= 0,.g/is 



where the time-dependent Fermi energy Ep{T) is deter- 
mined by 2 ^'^^'^Ulo p^^-^{T,Lo) = 1 when there is only 
one electron in the quantum dot. Due to the time period- 
icity introduced by the THz field, Ef{T) and M(r) are 
also periodic functions of T with period Tp. In Fig. 4, we 
plot the y-component of magnetic moment as a function 
of time. The magnetic moment is controlled by the RF 
strength, the frequency as well as the SOC and it can be 
qualitatively described by the factor A. The stronger the 
RF and the SOC are, the larger A and My become. On 
the other hand, the larger Q. is, the smaller A becomes 
and therefore the smaller My becomes. 
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FIG. 4: The average magnetic moment My versus the time for 
a single spin in QD at (a) E — 0.3 (solid curve), 0.5 (dotted 
curve), and 0.6 (dot-dashed curve) kV/cm for fixed SI = 4 
THz, and a_R = 3 x 10~® eV-cm; (b) an — 1 (solid curve), 
1.5 (dotted curve), and 3 (dot-dashed curve) xl0~^ eV-cm 
for fixed RF with = 4 THz and E = 0.6 kV/cm; (c) ^ = 4 
(solid curve), 4.6 (dotted curve), and 6 (dot-dashed curve) 
THz for fixed E = 0.6 kV/cm and Ofl = 3 x lO"'' eV-cm . 



In conclusion, we study the effects of the intense THz 
field on single electron spin in QD with the Rashba SOC. 
We calculate the single electron DOS at zero temperature 
and show that the DOS is greatly affected by the THz 
field and the SOC due to the dynamic Franz-Keldysh 
effect, the AC Stark effect and the side-band effect. It is 
shown that the joint effect of the THz field and the SOC 
can excite a THz magnetic moment which is controlled by 
the THz field strength, the THz frequency as well as the 
SOC. This provides a unique way to convert THz electric 
signals into THz magnetic ones which may be useful in 
full electrical magnetic resonance measurements. 
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